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THE  USE  OF  SOURCE-SINK  AND  DOUBLET  DISTRIBUTIONS 
EXTENDED  TO  THE  SOLUTION  OF  ARBITRARY  BOUNDARY 
VALUE  PROBLEMS  IN  SUPERSONIC  FLOW 
By  Max.  A.  Heaglet  and  Harvard  Lomax 

SUMMARY 

A  direct  analogy  is  established  between  the  use  of  source— sink 
and  doublet  distributions  in  the  solution  of  arbitrary  boundary  value 
problems  in  subsonic  wing  theory  and  the  corresponding  problems  in 
supersonic  theory.  The  concept  of  the  "finite  part"  of  an  integral 
is  introduced  and  used  in  the  calculation  of  the  improper  integrals 
associated  with  supersonic  doublet  distributions.  The  general 
equations  developed  are  shown  to  include  several  previously  pub¬ 
lished  results  and  particular  examples  are  given  for  the  loading  on 
rolling  and  pitching  triangular  wings  with  supersonic  leading  edges. 

INTRODUCTION 

The  problem  of  finding  pressure  distributions  over  airfoils  of 
arbitrary  shape  and  plan  form  or  of  finding  airfoils  which  have 
arbitrary  pressure  distributions  is  one  of  the  most  fundamental 
problems  in  aerodynamic  theory.  At  the  present  time  the  most 
important  and  satisfactory  approach  to  problems  of  this  type  is 
provided  by  the  methods  of  so-called  thin— airfoil  theory.  The 
essential  assumptions  in  this  theory  are  that  the  perturbation 
velocities  induced  by  the  airfoil  are  small  relative  to  the  free- 
stream  velocity  and  that  the  boundary  conditions  can  be  specified 
in  a  fixed  reference -plane . 

Under  the  assumptions  of  thin-airfoil  theory  the  theoretical 
analysis  of  a  problem  in  wing  theory  resolves  itself  into  the  task 
of  determining  the  solution  of  a  second— order  linear  partial 
differential  equation  with  prescribed  boundary  conditions.  In  the 
case  of  purely  subsonic  flow,  Laplace's  equation  in  three  dimensions 
must  be  considered,  while  in  purely  supersonic  flow  the  differential 
equation  which  arises  is  algebraically  equivalent  to  the  two- 
dimensional  wave  equation  of  mathematical  physics.  The  claseical 
solutions  of  these  two  equations  have  been  developed  along  two 
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distinct  lines:  first  by  use  of  orthogonal  functions  whioh  can  be 
derived  in  terms  of  the  boundary  conditions,  and  alternatively  by 
of  Green's  theorem  whioh  in  turn  utilizes  a  known  particular 
solution  of  the  partial  differential  equation  together  with  the 
given  boundary  conditions. 

One  particular  solution  associated  with  Laplace's  equation  and 
subsonic  aerodynamics  has  been  found  to  be  outstanding  in  its 
mathematical  usefulness  and,  when  identified  with  the  velooity 
potential,  has  a  physical  interpretation  which  can  supply,  in  direct 
application,  added  insight  into  the  nature  of  the  problem.  This 
function  is  referred  to  as  the  "fundamental  solution"  and  can  be 
developed  from  the  concept  of  a  so— called  souroe.  A  concomitant 
development  to  the  souroe  potential  is  the  doublet  potential,  and 
appropriate  distributions  of  these  functions  are  known  to  be 
auff ic lent  for  the  solution  of  all  problems  in  subsonic  wing  theory.' 

The  extension  of  the  use  of  the  fundamental  solution  to  problems 
in  purely  supersonic  flow  introduces  mathematical  difficulties  which 
differ  essentially  from  those  encountered  at  low  speeds.  Both  the 
source  and  the  doublet  potentials  possess  singularities  on  their 
conical  oharaoteristio  surfaoes  or  Mach  cones  and,  in  the  case  of 
the  doublet,  the  singularity  is  of  higier  order  than  oan  be  treated 
by  elementary  mathematical  methods.  In  the  historical  development 
of  the  solutions  of  the  wave  equation  this  trouble  was  oircumvented 
by  replacing  the  souroe  potential  by  other  particular  solutions  of 
the  differential  equation.  As  an  example,  Yolterra  (reference  1) 
introduced  the  integral  of  the  fundamental  solution  and  in  that  way- 
reduced  the  order  of  the  singularities  involved.  The  analytical 
development  of  Yolterra*  s  theory  presents  no  inherent  difficulties 
(e.g.,  reference  2)  but  the  physical  significance  of  the  particular 
solution  is  lost,  the  direot  analogy  with  subsonio  theory  no  longer 
exists,  and  a  certain  amount  of  mathematical  inefficiency  arises 
since,  after  using  the  integral  of  the  souroe  potential,  it  is 
found  neoessary  to  resort  at  the  end  of  the  analysis  to  taking  a 
final  derivative. 

In  this  report,  following  methods  introduced  by  Hadamard 
(reference  3),  a  general  solution  to  the  thin-airfoil  problem 
in  superscmio  theory  will  be  given  in  terms  of  the  distribution  of 
sources  and  doublets  over  the  given  reference  plane.  Furthermore, 
a  discussion  of  the  nature  of  the  boundary  values  required  will  be 
given*  For  properly  set  problems  in  wave  theory  it  has  been  found 
necessary  to  specify,  usually,  both  the  required  function  and  its 
derivative  with  respect  to  time  along  the  boundary  considered.  In 
aerodynamio  applications  of  the  wave  equation  associated  with 
lifting  surface  theory  and  thickness  distributions  it  will  be  shown 
that  only  a  knowledge  of  the  unknown  function  or  its  normal 
derivative  along  the  boundary  is  needed  since  a  relationship  between 
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the  two  functions  will  be  established  cm  the  boundary  surface . 

In  the  theoretical  portion  of  the  report  a  brief  presentation 
will  be  made  of  the  differential  equations  involved  and  the  two  forms 
of  the  fundamental  solution.  An  outline  is  then  given  of  the  types 
of  boundary  value  problems  encountered  and,  since  the  purpose  of  the 
report  is  to  extend  the  concepts  of  thin-airfoil  theory  which  are 
used  in  subsonic  theory  to  problems  arising  in  supersonic  theory,  a 
discussion  will  be  given  of  the  subsonic  development  as  a  basis  for 
the  analogy  which  exists  between  the  methods  of  solution  corresponding 
to  the  two  regimes  of  flow.  In  the  discussion  of  the  purely  supersonid 
case  it  will  be  shown  that  the  introduction  of  the  concept  of 
"finite  part"  will  provide  a  technique  whereby  the  Improper  integrals 
arising  from  the  use  of  doublets  may  be  evaluated  in  a  straight¬ 
forward  manner.  The  applications  of  the  theoretical  developments 
will  include  the  rederivation  of  same  previously  published  results 
and  will  also  contain  the  calculation  of  load  distributions  for 
rolling  and  pitching  triangular  wings  with  leading  edges  swept  ahead 
of  the  Mach  cone  from  the  vertex  of  the  triangles. 
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B  arbitrary  region  of  Integration 

S  surface  enclosing  region  B 

u,v,w  perturbation  velocities  in  direction  of  X—  ,  T—>  and  Z— axes, 
respectively 

T  f re e-stream  velocity 

I,T,Z  Cartesian  coordinates  In  original  spaoe  variables 
XfjfZ  transformed  system  of  coordinates 

0  J  |  M2-lj 

c  Infinitesimal  used  In  analysis 

X  surface  along  which  Btream  enters  induoed  field  of  wing 

t)  conormal  to  arbitrary  surfaoe 


n1#n2,oa  direction  cosines  of  conormal 

a  variable  representing  either  acceleration  potential, 

Telocity  potential,  or  any  of  the  three  perturbation 

Telocity  components 

T  surface  on  whioh  boundary  conditions  are  given 

$  perturbation  velocity  potential 

A  variable  representing  either  acceleration  potential, 

velocity  potential,  or  any  of  the  three  perturbation 

velocity  components 
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differential  operator 
sign  denoting  "finite  part" 


of  integral 


Subscripts 

U  SU^glpt  den°tinS  Value  of  liable  on  upper  surface  of 

1  8U™  denotinS  Talue  of  variable  on  lower  surface  of 

1  subscript  denoting  variable  of  integration 

"  subscript  on  r  denoting  fundamental  solution  in  supersonic 


Superscript 

superscript  denoting  value  of  variable  on  opposite  side 
of  t  from  filed  point  (x,y,z) 


THEORETICAL  DEVELOPMENT 

Linearized  Equations  and  Boundary  Conditions 

The  linearization  of  the  second-order  differential  equation  for 
compressible  fluid  flow  is  developed  under  the  assumptions  of  thin- 
alrf oil  or  small-perturbation  .theory.  If  the  velocity  vector  of  the 
free  stream  Is  parallel  to  and  in  the  direction  of  the  positive 

X-axis,  the  resulting  differential  equation  Is  expressible  in  the 
form 


<»«■>  g  ♦ 

ax2 


+ 

dY2  +  dZ2 


(1) 


where  fl  represents  a  velocity  potential,  acceleration  potential,  or 
any  one  of  the  perturbation  velocities  while  M  is  the  constant 
value  of  the  free— stream  Mach  number.  Assuming  the  plane  of  symmetry 
of  the  airfoil  to  lie  in  the  XY  plane,  the  boundary  conditions 
associated  with  equation  (l)  are  given  for  Z=0.  Moreover,  if 
u,  v,  and  w  represent,  respectively,  the  perturbation  velocity 
components  along  the  X,  Y,  and  Z-axes,  and  if  the  velocity  of  the 
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free  stream  is  V,  the  direction  cosines  of  any  stream  line  are 
proportional  to  the  point  functions  V  +  u(X*  Y,  Z),  v(X,  Y,  Z) , 

and  w(X,  Y,  Z)  while  pressure  coefficient  —  is  given  by  the 
relation 

Zip  2u  , 

T  =  “  “  1 

Detailed  discussions  of  these  results  may  be  found  in  reference  4. 
Introducing  the  affine  transformations 


1  =  X 

y  =  «/±(l-M2)  Y 

z  =  Z 

where  the  signs  under  the  radicals  are  chosen  so  that  real  values 
result,  it  follows  that  in  the  subsonic  case  (M<l)  equation  (1) 
reduces  to 

ifs  ♦  =  0  < 

dx2  dy2  dz2 

while  the  supersonic  case  (M>l)  yields 

a2n  a2^  _  S2^  _  ^ 

dx2  dy2  dz2 

The  fundamental  solution  associated  with  equation  (3)  i0 

i  =  [(x-xj.)2  +  (y-yi)2  +  (z-zi)  1  2  ( 


or,  in  terms  of  the  original  space  variables. 


=  [(X-Xx)2  +  P2(Y-Yl)2+  p2(Z-Zi)£] 


where 


p2  =a-M2) 
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When  the  wave  equation  is  to  he  considered  the  fundamental  solution 
takes  the  form 

—  =  [(x-sa)1 2  -  (y-yi)2  -  (z-zi)2]  2 

rc 
or 

_  i 

—  =  [  (X-X!)2  -  02(Y— Yi) 2  -  02(Z-Z!)2]  2 

rc 

where 

02  =  (M2  -  1) 

These  fundamental  solutions  represent,  respectively,  in  subsonic 
and  supersonic  flow  the  velocity  potentials  at  the  point  (x,y, z)  or 
(X,Y,Z)  of  unit  sources  situated  at  the  point  C*i,yi,zi)  or 
(xi^Y1,Z1).  The  velocity  potential  of  a  doublet  may  be  obtained  by 
taking  a  directional  derivative  of  the  source  potential,  the 
direction  of  the  axis  of  the  doublet  coinciding  with  the  direction 
along  which  the  derivative  is  taken.  These  two  functions  will  be 
seen  to  be  of  paramount  importance  when  Green* s  theorem  is  applied 
to  the  given  boundary  conditions. 

It  remains  now  to  mention  the  types  of  boundary  conditions  which 
appear  in  problems  associated  with  wing  theory.  As  a  convenience  to 
the  development  of  the  theory  the  normalized  forms  (equations  (3) 
and  (4))  of  equation  (l)  will  be  used  and  boundary  conditions  will  be 
assumed  known  with  respect  to  the  x,y,z  coordinate  system.  Retransfor— 
mation  totheX,Y,Z  system  of  axes  can  be  made  quite  Bimply  wherever 
needed  in  application.  In  order  to  define  the  boundary  conditions, 
two  subscripts  will  be  introduced:  the  first,  u,  denotes  the  value 
of  the  required  function  on  the  upper  surface,  that  is  the.  limit  of 
the  function  as  z  approaches  zero  from  the  positive  direction;  the 
second,  I,  denotes  the  value  on  the  lower  surface,  that  is,  the 
limit  of  the  required  function  as  z  approaches  zero  from  the  negative 
direction. 

Using  these  definitions  the  three  boundary  value  problems  of 
principal  interest  can  be  defined  as  follows: 

1.  Symmetrical  nonlifting  airfoil.—  In  this  case  wu  =  wj  =  0 
over  all  of  the  xy— plane  exoept  for  the  region  occupied  by  the 
airfoil  where  2wu  =  — 2w^  =  Aw  =  f(x,y)  the  function  being  determined 

by  the  geometry  of  the  wing.  Over  all  of  the  xy— plane,  Au  =  0. 

2.  Lifting  plate  with  specified  loading.—  It  is  given  that 
Au-=  uu  —  ux  =  0  over  the  xy-plane  except  for  the  region  occupied 


(6) 

(6a) 
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by  the  airfoil  where  Au  =  f(x,y),  the  function  being  determined  by 
the  specified  loading.  Moreover,  Aw  =  0  everywhere. 

3.  Lifting  plate  with  specified  camber,  twist,  and  angle  of 
incidence.—  Over  the  xy— plane  Aw  =  0  everywhere.  And,  except  for 
the  region  occupied  by  the  airfoil,  Au  =  0.  Over  the  region  occupied 
by  the  airfoil  w  =  f (x,y)  where  f (x,y)  is  determined  by  the  given 
camber,  twist,  and  angle  of  incidence. 

It  should  be  pointed  out  that  the  problems  considered  here  differ 
from  the  usual  type  of  boundary  value  problem  encountered.  In  the 
so-called  Dirichlet  or  Neumann  problems,  which  arise  in  connection 
with  Laplace fs  equation,  the  value  of  the  normal  derivative  of  the 
function  or  of  the  function  itself  is  specified  along  the  boundary 
while  the  Cauchy  problem  for  second— order  partial  differential 
equations  involves  the  knowledge  of  both  the  function  and  a  derivative. 
Except  for  one  case  in  the  aerodynamic  problems  listed  above,  no 
absolute  values  are  given  but  rather  the  Jump  in  the  value  of  the 
function  along  the  boundary  is  prescribed. 


Boundary  Yalue  Problems  in  Purely  Subsonic  Flow 

Since  the  purpose  of  this  report  is  to  extend  the  concepts  of 
thin-airfoil  theory  which  are  used  in  subsonic  theory  to  problems 
arising  in  supersonic  theory,  some  discussion  of  the  former  will  be 
given  to  provide  lucidity  as  well  as  to  furnish  a  basis  for  the 
analogy  which  will  be  shown  to  exist  between  the  methods  of  solution 
arising  in  the  two  regimes  of  flow. 

The  method  whereby  the  solutions  of  the  given  problems  can  be 
effected  is  provided  by  Green's  theorem  which  relates  a  volume 
integral  over  a  region  R  to  a  surface  integral  over  the  surface 
S  enclosing  R.  If  a,  ft  are  any  two  functions  which,  together 
with  their  first  and  second  derivatives,  are  finite  and  single— valued 
throughout  R,  then  for  the  subsonic  case 


n  ^ = f£  f  *■ a  v2n  ~ n  ^ 


where  the  Laplacian  operator. 


V2 


^  .  d2  .  d2 
=ZT*+  'Sp  +  5SS' 


is 


(7) 


introduced  and  the  directional  derivatives  on  the  left  side  are 
taken  along  the  normal  n,  drawn  inward,  to  the  surface  S.  ^ 
Identifying  now  the  function  a  with  the  fundamental  solution  — 

and  specifying  that  ft  satisfies  Laplace's  equation,  equation  (7) 
simplifies  to  give 
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where 

i  =[  (x-ii)2  +  (y-yi)2  +  (z-Zi)2r  * 

The  variables  of  integration  in  the  equation  are  x*,  yx,  z±}  while 
7 ,  z  are  the  coordinates  of  a  point  P  either  inBide  or  outside  of 
the  region  of  integration. 

If  the  point  P  is  assumed  to  be  inside  the  region  of  integration, 

it  is  evident  that  the  function  i  becomes  infinite  at  P.  and  it  is 

r 

necessary  to  exclude  this  point  from  the  region  if  formula  (8)  is  to 
apply.  Describing  a  spherical  surface  E  with  radius  e  about  the 
point  P,  and  considering  the  integral,  over  the  two  surfaces  E  and  8 
which  enclose  the  region,  it  can  be  shown  that  in  the  limit  as 
£  — •>  0  equation  (8)  becomes 


(x,y,z) 


(9) 


The  physical  significance  of  this  last  relation  follows  Immediately: 

the  term  —  represents  a  fluid  source  and  the  term 

r  dn 

represents  a  doublet  with  its  axis  lying  along  the  normal  to  S,  both 
source  and  doublet  being  situated  at  the  surfaoe  point  Xj.,  yi,  zi. 

The  value  of  the  function  H  at  the  point  x,  y,  z  is  therefore  given 
as  an  integral  of  a  source  and  doublet  distribution,  the  strengths  of 
the  two  being  determined  directly  from  the  respective  boundary  values 

of  fl  and  ££ 
dn  ' 

Equation  (9)  expresses  the  value  of  ft  in  terms  of  the  surface 

values  of  ft  and  hut  this  relation  does  not  Imply  that  a 

on 

knowledge  of  both  these  variables  is  necessary  for  the  determination 
of  ft  .  As  can  be  shown  easily,  another  condition  may  be  established 
which  relates  the  two  surface  values.  Applying  equation  (8)  to  the 
case  where  P  lies  outside  the  region  of  integration,  it  follows  that 


the  integral  is  equal  to  zero  and  that  ft  and 
are  functionally  dependent. 


on  the  surface 
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Sufficient  information  is  now  at  hand  to  provide  a  solution  for 
the  thin-airfoil-houndary  value  problems.  Consider  the  region  R 
hounded  by  the  xy— plane  and  a  hemispherical  dome  of  infinite  radius 
lying  above  this  plane.  For  all  problems  to  which  the  results  will 
be  applied,  the  value  of  ft  may  be  assumed  equal  to  zero  at  infinity. 
The  contribution  of  the  surface  integral  over  the  hemisphere  is  thus 
zero  and,  from  equation  (9) 


dJCi  dyx 


where  the  integration  extends  over  the  entire  plane.  The  directional 
derivatives  are  necessarily  in  the  direction  of  the  positive  Zr-axis 
and  subscripts  are  again  introduced  to  denote  conditions  on  the  upper 
side  of  the  plane.  Keeping  P  fixed  and  integrating  over  the  lower 
side  of  the  xy— plane ,  it  follows  that 

where  the  negative  direction  of  the  normal  may  he  ignored  since  the 
integral  is  equal  to  zero.  Subtracting  these  two  equations  gives 
the  expression 

“  (w>  -  -  ^  (10> 

the  integral  extending  now  only  over  the  area  r  for  which  the 
integrand  does  not  vanish.  Equation  (10)  is  the  basic  equation  from 
which  all  solutions  in  subsonic  wing  theory  will  be  developed.  It 
should  be  pointed  out  that  the  derivation  proceeded  from  the 
assumption  that  the  point  (x,y,z)  lay  above  the  xy-plane.  When 
(x,y,z)  lies  below  the  xy— plane ,  however,  the  derivation  can  be 
carried  through  in  exactly  the  same  manner.  Such  a  development 
reveals  that  equation  (10)  is  general  and  that  no  restriction  need 
be  imposed  on  the  position  of  (x,  y,  z)  relative  to  the  reference 
plane. 

As  a  particular  application  of  equation  (10)  consider  a  thin 
symmetrical  airfoil  at  zero  angle  of  attack  and  set  ft  =  ®  where 
$  is  the  perturbation  velocity  potential.  Conditions  of  symmetry 

demand  that  0^  =  =  ft}  while  =  wu  and  =  wj. 
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Thus,  if  wu  —  wj  =  Aw 

*  ■  “  hfl*  ; 311 471  (11> 

and  the  velocity  potential  is  given  by  a  distrihution  of  source 
potentials.  This  distrihution  can  he  immediately  related  to  the 
slope  of  the  basic  section  hy  means  of  the  equation 

2!h  -  _  IL  -  f^C) 

V  ~  ~  V  \dx/u 

The  symmetric  airfoil  can  also  he  treated  hy  replacing  A  hy 
the  perturbation*  velooity  w  and  in  the  case  of  the  thin  lifting 
surface  with  given  loading  the  function  A  can  he  set  equal  to  u. 
Employing,  respectively,  conditions  of  symmetry  and  irrotationality, 

it  follows  that  ^2^—  vanishes  and,  after  setting  A  A=  Au  —  A 
dz  dz 

equation  (10)  becomes 

A  (x,y,z)  *  i-  yyk  a|-  0  dXi  dyj.  (12) 


Boundary  Value  Problems  in  Purely  Supersonio  Plow 

Applications  of  Green’s  Theorem.—  The  problem  to  he  discussed 
at  this  point  is  the  extent  to  which  an  analogue  to  equation  (10)  can 
he  developed  for  supersonic  flow  fields.  The  first  step  in  the 
presentation  is,  once  more,  the  introduction  of  Green's  theorem  for 
equation  (l)  after  it  has  been  modified  to  the  form  given  hy  equa¬ 
tion  (4).  Rnploying  the  operator 

n2  =  £L_a2  -*1 

dx2  cly2  dz2 


Green's  theorem  now  becomes 


(13) 


% 
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where  V)  is  the  so-called  conormal  to  the  surface  S  and  has 
direction  cosines  equal  to  V>2,  t>3  such  that 

=  ”2-  "s  =  ”3  (14) 

where  ni>  n.2,  113  are  the  direction  cosines  of  the  normal  to  the 

surface  S  (fig.  l) .  (The  conormal  at  any  point  xlf  yi,  Zi  of  a 
surface  is  the  mirror  image  in  the  plane  x  =  x*  of  the  normal 
through  the  same  point.)  If  a  and  n  are  perfectly  arbitrary 
functions,  aside  from  satisfying  the  usual  conditions  of  single— 
valuedness,  etc.,  equation  (13)  represents  an  identity  and  this,  fact 
will  be  useful  at  a  later  time.  For  -Immediate  purposes,  however,  a 
and  Q  will  be  chosen  as  solutions  of  the  differential  equation  under 
consideration  so  that 

Ua  =  =  0 

and,  consequently. 


a  —  -  0  — ^)  ds  =  0  (15) 

So  bxy 

The  use  of  equations  (13)  and  (15)  depends  upon  an  understanding 
of  the  physical  nature  of  supersonic  flow  fields.  The  essential 
feature  of  such  flow  is  the  presence  of  Mach  cones  which  correspond 
to  the  characteristic  cones  arising  in  the  mathematical  study  of  the 
wave  equation.  In  accordance  with  these  concepts  a  disturbance  in 
the  flow  field  can  affect  the  flow  only  within  its  aftercone,  that 
is,  the  cone  with  vertex  at  the  point  of  disturbance  and  with  axis 
extending  in  the  direction  of  the  undisturbed  stream  velocity 
vector;  conversely,  a  point  in  the  flow  field  can  be  affected  only 
by  disturbances  which  emanate  from  points  within  its  forecone. 

When  the  disturbances  are  generated  by  a  wing  it  is,  moreover, 
necessary  to  speak  more  specifically  about  the  nature  of  the  leading 
edge  of  the  wing.  For  all  cases  considered  here  the  assumption  will 
be  made  that  the  plan  form  is  a  polygon,  that  is,  is  composed  of 
straight  line  segnents.  If  the  wing  is  swept  ahead  of  the  foremost 
Mach  cone,  the  cones  arising  at  the  leading  edge  will  have  as 
envelope  a  wedge— shaped  surface  passing  through  and  extending  back 
from  the  leading  edge,  while  if  the  wing  is  swept  back  of  the 
foremost  Mach  cone  this  cone  will  be  the  surface  along  which  the 
air  first  experiences  perturbations  or  disturbances.  Thus,  a 
point  P  with  coordinates  x,  y,  z  is  affected  by  all  disturbances 
lying  within  its  forecone  T  and  at  the  same  time  behind  the  forward 
surface  X,  the  nature  of  the  latter  surface  being  dictated  by  the 
leading  edge.  In  figures  2(a)  and  2(b)  these  surfaces,  along  with 
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the  disturbance  plane  t,  are  indicated  for  two  different  wing  plan 
forms.  In  the  applications  of  equation  (13)  the  volume  integral  is 
limited  to  the  portion  of  space  common  to  the  surfaces  r,  A,  and  t 
and  the  surface  integral  involves  a  discussion  of  conditions  on  these 
surfaces. 

Up  to  this  point  the  analogy  between  the  subsonic  and  supersonic 
cases,  insofar  as  the  use  of  Green’s  theorem  is  concerned,  is  quite 
apparent.  The  principal  difference  which  occurred  was  brought  about  by 
the  use  of  the  true  normal  in  the  subsonic  field  together  with  the  fact 
that  the  xy-plane  was  covered  by  a  hemispherical  dome  of  infinite  radius 
whereas.  In  the  supersonic  field,  the  concept  of  the  conormal  was 
introduced  and  the  volume  to  be  considered  was  that  enclosed  within 
^  finite  region.  In  continuing  the  analogy,  however,  far  more 
formidable  obstacles  arise.  To  begin  with,  the  discussion  of  a  and 
fi  over  the  surface  in  the  subsonic  case  was  relatively  simple.  Thus, 
with  no  limitations  of  generality  G  could  be  assumed  zero  at 
infinity  and  AG  was  specified  completely  in  the  xy— plane.  But  in 
the  supersonic  case,  although  AG  can  be  assumed  known  In  the  xy— plane 
and,  as  will  be  seen  later,  G  may  be  evaluated  on  the  forward  boundary 
of  the  region,  nothing  is  known  of  G  on  the  forecone  T.  Hence  a 
must  be  chosen  properly  so  that  the  knowledge  of  G  is  unnecessary  on 
r.  The  most  obvious  choice  of  <j  would  be  a  particular  solution  of 
equation  (4)  which  would  make  a  =  0  on  T  and  this  is  in  fact  the 
choice  used  by  Yolterra  (reference  l)  and  applied  to  aerodynamic 
problems  in  reference  2.  However,  if  the  analogy  is  to  be  maintained 
the  choice  of  a  is  not  arbitrary  but  must  be  the  three-dimensional 
supersonic  source  corresponding  to  the  fundamental  solution 

-  in  subsonic  theory.  But  such  a  solution,  —  = 
r  rc 

[  (x-xx)2  -  (y— yi)2  -  (z-zx)2]-  2  becomes  infinite  along  the 

forecone  T  which  has  the  equation  (x— X1)2  —  (y— yi)2  —  ( z— z x ) 2  =  0; 

It  is  Just  thiB  difficulty  which  apparently  invalidates  any  further¬ 
ance  of  the  analogy  and  the  prediction  in  advance  of  an  aerodynamic 
shape  from  a  distribution  of  sources  and  doublets  in  supersonic  flow. 
However,  it  is  also  precisely  this  difficulty  which  is  overcome  by 
Hadamard’s  general  methods. 


Extension  of  analogy  by  Hadamard’s  Method.-  The  full  develop¬ 
ment  of  Hadamard’s  methods  cannot  be  given  here,  but  a  rough  sketch 
of  his  reasoning  will  perhaps  be  useful.  The  basis  of  his  arguments 
stems  from  equation  (13).  First  it  is  admitted  that  the  right-liand 
side  of  equation  (13)  will  tend  to  infinity  as  the  surface  S 
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approaches  r  so  that  —  is  not  a  regular  solution  to  1  1  Q  =  0 

^■*c 

on  T.  However,  as  has  teen  mentioned,  equation  (13 )  still  must  hold 
whether  or  not  a  or  ft  satisfy  the  wave  equation  and  thus  it  still 
provides  an  equality.  Hence,  'if  the  surface  integral  tends  to 
infinity  so  also  must  the  volume  integral.  Further,  equation  (13) 
implies  that  these  infinite  portions  Just  cancel  since  the  difference 
of  the  two  integrations  must  always  give  zero.  To  deal  with  such  a 
problem  quantitatively  by  the  usual  mathematical  techniques  would 
require  the  study  of  a  limiting  process  for  each  new  boundary  value 
problem.  Hadamard's  contribution  was  the  introduction  and  Justifica¬ 
tion  of  a  concept  which  removed  the  necessity  for  studying  the 
infinite  portions  involved.  This  concept  is  best  presented  by 

means  of  a  new  notation,  thus  the  sign  |  is  used  and  is 

to  be  read  "the  finite  part  of." 

Using  this  concept  it  is  possible  to  show  that  if  a  were  set 
equal  to  — ,  then  equation  0-3)  could  be  written 


so  that  the  "finite  parts"  of  eaoh  side  of  the  equation  would  be 
equal.  Such  a  notation  would,  of  course,  in  general  be  meaningless 
since  in  discarding  arbitrarily  a  part  which  tended  to  infinity  it 
would  be  possible,  by  proper  combinations,  to  obtain  as  a  remainder 
any  finite  value.  The  fact  is,  however,  that  the  integrals  involved 
in  equation  (l6)  tend  to  infinity  only  at  a  limit  of  the  integration 
and  this  limit  always  involves  the  forecone  T .  It  was  consequently 
possible  to  devise  a  manipulative  technique  to  handle  equation  (l6) 
without  considering  the  singularities  individually.  It  might  be 
mentioned,  without  stressing  the  correspondence,  that  a  treatment  of 
improper  integrals  is  also  employed  in  the  use  of  Cauchy's  principal 
value.  In  subsonic  thin— airfoil  theory  and  lifting— line  theory 
integrals  of  the  latter  type  are  well  known  in  the  form 
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II  = 


f(x0)dx0 

I-Xp 


o<  x  <c 


Ii  certainly  tends  to  infinity  as  x0  approaches  x  hut  the  use 
of  Cauchy's  principal  value  allows  the  very  large  values  of  the 
integrand  obtained  when  x0  is  on  either  side  of  x  to  Just  cancel 
in  such  a  way  that  Ix  is  finite  and  unique.  So  again  the 
integral 


la  = 


xo  A(x)dx 


(xQ-x) 


3/2 


is  finite  and  unique  and  given  by  Hadamard  in  the  form 

xo 


-f. 


A(x)-A(x0) 
(xp-x) 


3/2 


dx 


2A(x0) 

(xp-a) 


1/2 


It  is  actually  possible  to  generalize  the  idea  of  "finite  part"  to  the 
case  when  the  exponent  in  the  denominator  is  of  the  form  j  +  —  where 
j  is  a  positive  integer  hut  such  a  generalization  is  not  needed  for 
aerodynamic  applications  and  will  therefore  he  omitted. 

In  actual  calculation*  the  evaluation  of  the  integral  I2  can 
he  shortened  considerably.  Thus*  if  the  indefinite  integral  of 


/ 


A(x) dx 

U O^F572 


is  written  in  the  form  f(x)-+C  then 


I2  = 


11m 

x  — >  x, 


X 


1  A(xi)-A(x0) 
(xn-xx)3/2 


dxi  — 


(xQ-a)l/= 


] 


lim  F(x)+C-F(a)-C— 

:  — >  x0  L 


2A(x0) 

(xo-x) 


1/2 


It  follows  that  if  C  Is  chosen  so  that 

2A(x0) 


C  =  lim 
x  — >  x, 


y(x) 
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then  the  egression  for  I2  may  be  written  I2  =  -[P(a)+C]  . 

When  C  is  chosen  in  this  manner,  the  notation  for  the  calculation 
may  thus  be  modified  to  the  form 


p*Q  A(x)dx  _  p*  A(x)cb 

i.  (x0-x)3/e“Ja  (Io_z)3 


A(x)dx 


where  the  asterisk  indicates  that  the  upper  limit  is  not  substituted 
into  the  indefinite  integral,  F(x)  +  C, 

The  technique  for  the  calculation  of  the  finite  part  has 
therefore  been  reduced  to  three  simple  steps:  first,  the  indefinite 
integral  r(  x)  +  C  is  determined,  second,  the  constant  C  in  tjie 
indefinite  integral  is  evaluated  hy  means  of  a  limiting  process, 
third,  the  lower  limit  of  the  integral  is  substituted  into  the 
indefinite  integral  and  a  minus  sign  prefixed.  As  an  example, 
consider  the  integral 


'a  (x02-x2)3/2  (x^-x2)3/2 


In  this  case  F(x)  +  C  =  — — i— — —  +  0  and 

(x0^-x2) 


1  1 

10  (2Xo)' v*  (xo-x)1^  (xo^x2)  V* 


so  that,  finally. 


(xoa-x2)3/2 


=  -EF(a)  +  C]  = 


(Xf,3-®-2)  V2 


With  the  aid  of  this  artif ioe  the  analogy  between  the  subsonic 
and  supersonic  cases  can  be  continued  with  relative  ease.  Thus,  in 
equation  (16)  the  rightr-hand  member  is  zero  provided  we  exclude  the 
point  P  from  the  volume  of  integration.  This  can  be  done  most  easily 
by  limiting  the  integration  to  the  ii  =  constant  plane,  a  distance  e 
upstream  from  P,  The  portion  of  this  plane  intersected  by  the  cone, 
and  thus  the  section  over  which  the  integration  must  be  carried,  will 
be  denoted  by  z  (fig.  3). 
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As  drawn,  figure  J  shows  a  cross  section  in  a  Jx  -  constant 
plane  for  the  special  case  when  P  is  located,  directly  behind  and 
above  the  foremost  disturbance.  Applying  equation  (l6)  to  the 
regions  above  and  below  the  disturbance  surface  r  (plane  of  the 
airfoil)  yields  the  two  equations 


If, 


A+r+T 


0 


d  /  l\  1_  dft  " 
do  \rc/  rc  do  - 


dS 


(17) 


and 


IL 


A+r+T 


L  1 

'V) 

_  _1_ 

dft»  ’ 

do‘ 

rc 

do'  - 

dS  =  0 


(18) 


where  the  prime  indicates  the  surface  value  of  ft  on  the  opposite 
side  of  t  from  P. 

The  integration  over  E  can  be  computed  for  e  very  small. 

For  convenience,  consider  P  to  be  the  origin;  then  it  follows  that 
since  the  conormal  is  in  the  x*  direction  and  the  area  element  can 

be  written  7  <Ly  d0  where  6  =  arc  tan  —  and  7  =  Jjx2  +  Zj.2,  the 

y  1 

right  side  of  equation  (17)  will  give 


=  11m  ■^7rfl(x,j ,z)  f  g  +27 r  f  - Z_J7_  |  =  airn(x.y.*) 

e_>0L  4>  Ce  _  7  ] '  bxx  J0  /eg'Z~ yg~  J 

Hence  the  value  of  Q  at  the  point  P,  Q  (x,y,z)  can  he  determined 
from  equation  (17)  vith  the  restriction  implied  by  equation  (l8). 

Further,  since  only  the  "finite  part"  is  considered,  the  integration 
over  T  yields  zero  and  the  two  equations  combine  to  give 


ft  (x,y,z)  = 


1_ 

2k 


rr 

1  /an  afi,v\  a  / 

J  JT 

-  rc  \  dz  bz'  dz 

1_ 

2k 


fl 


a 

,a  — 

L  So  \r. 


1  3fll  l 

- as - 

rc  do  J  2k 


IL 


Q 


a 

a-o* 


aa 

au*  j 

(19) 


dS 


The  only  remaining  difference  between  the  subsonic  solution  for  the 
distribution  of  sources  and  doublets,  equation  (10),  and  the  super¬ 
sonic  solution,  equation  (19),  is  the  integration  over  surface  X. 
The  detailed  discussion  of  the  contribution  of  the  surface  integrals 
over  X  will  be  deferred  to  Appendix  A.  It  must  suffice  for  the 
present  to  remark  that  in  all  applications  the  integrals  over  X  in 
equation  (19)  are  either  zero  or  combine  to  give  zero.  It  therefore 
follows  that 
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fi(x,y,z) 


1_ 

2* 


oQu  __ 
d  z  bz 


bz 


0 


(20) 

and  the  complete  analogue  to  equation  (10)  has  been  established  through 
the  use  0/  the  concept  of  finite  part. 


APPLICATIONS 

Interpretation  of  Previous  Results 

As  a  means  of  indicating  the  various  problems  to  whioh  equation  (20) 
can  be  applied,  three  previously  published  results  will  be  discussed. 
These  applications  include,  first,  the  expression  for  the  perturbation 
velocity  potential  of  a  symmetrical  non lifting  airfoil  (reference  6), 
second,  the  calculation  of  pressure  distribution  over  a  semi— infinite 
wedge  with  leading  edge  swept  back  of  the  foremost  Mach  cone 
(reference  5)  and,  third,  the  integral  equation  method  for  determining 
the  load  distribution  over  a  lifting  surface  of  arbitrary  slope 
(reference  7) • 

As  in  the  case  of  equation  (11)  for  subsonic  flow,  let  ft 
represent  velocity  potential  9  and  consider  the  case  of  a  symmetrical 
wing  at  zero  angle  of  attack. 


Then  — — ^  =  wu  and  — 
dz  dz 


=  Wj 


where  wu  and  w^ 


are 


induced  vertical  velocities  on  the  upper  and  lower  aurfaceSj 


respectively.  Moreover,  $  —  <D  =  0  for  the  symmetrical  case  so 

u  if 

that,  since  wu  —  Wj  =  2  wu. 


9 


dyx 


The  integral  in  this  equation  is  finite  at  T  so  the  finite  part 
sign  may  be  disregarded  and 


9 


dxx 


dyx 


(21) 
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This  equation  agrees  with  results  given  by  Puckett  in  reference  6. 

As  another  example  consider  the  solution  used  h y  R.T.  Jones 
in  reference  5  for  a  nonlifting  symmetrical  wing.  Setting  Q 
equal  to  w  in  equation  (20)  and  using  the  fact  that  w  and  u 
are  related  by  the  expression 


u 


w  dz 


it  follows  that 


u 


dxi  dyi 


(22) 


For  a  wedge  swept  behind  the  forward  Mach  cone  and  having  as  the 
equation  of  its  leading  edge  the  relation  Xi  =  myi,  the 
expression  for  u  may  be  written  in  the  form 


u 


it  <5x 


z  dxi 

[(x-Xi)2  -  (y-yx)2  -  z2]3/2 


where 


Performing  the  integration  with  respect  to  xlt  it  follows  that 


u 


=  Si  L 

«  dx 


zdz 


U  -  £0  dyi 


((7-yi)2  +  Z2]  J  (x  -  |i)2-  (y— yi)2  -  z2 


and,  after  reversing  the  order  of  integration 
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w„  d  r>Yl  n 


_ -  (x  -  tt)  zdz 

C  (y-yi  )s+zs]  J(x  -  ?±)2- 


-)2-(y-7i)2-z2 


=  k 


fui  rY_^ 

2rt  dx  J0 


In 


\/  (x  ~  “  (y-7i)2  -  z2  -  k  (x  - 

L  1  U  -  |i)2  -  (y-yi)2  -  z2  +  k  (x  -  Ji)  -J 


dyi 


where  k  =  1  for  x>—  and  -1  for  x<—  .  Taking  the  partial 

derivative  with  respect  to  x  and  noting  that  the  value  of  the  logarithm 
at  the  upper  limit  is  zero,  the  value  of  the  induced  velocity  is 


and  integration  yields  the  final  result 


wu  m  (x  -  my)  -J  (1  -  m2)(x2  -  y2  -  Z2) 

In - j—  ■  —  - . . .  ■  - - 

*/l  —  m2  J (mx  -  y)2  +  z2  (1  —  m2) 


Denoting  pressure  coefficient  —  by - and  setting  —  equal 

q  7  7 

to  l  — 


£p  _  2  ( iz 

q  rt 


the  slope  of  the  surface,  this  may  be  written 
(x  -  my)  -  J  {\  -  m2)  (x2  -  y2  - 


In 


z2) 


"(mx  -  y)2  +  z2  (1  -  m2) 


(23) 


Equation  (23)  gives  the  pressure  coefficient  at  any  point  in  the 
field  produced  by  a  wedge  swept  behind  the  Mach  cone.  When  z  is 
set  equal  to  zero  the  pressure  distribution  over  the  wedge  itself  is 
determined  and  the  equation  corresponds  exactly  with  equation  (12) 
of  reference  5* 
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When  loading  la  to  he  prescribed  over  a  thin  lifting  surface,  fl 
may  be  assumed  equal  to  the  perturbation  velocity  u.  A  direct 
consequence  of  this  assumption  is  that  in  equation  (20) 

Sz  dz 


since,  from  conditions  of  irrotationality, 

dz  dx  dx  dz 

By  definition 

.  2uI 

q  V  '  q  V 

and  load  distribution  in  coefficient  form,  ACp,  is  given  by  the 
relation 


ap  .£l_£Si 

9  9 


so  that 


2u,  2u  An 

A3p=--Yi  +  ^  =  -2^r 


Equation  (20)  can  therefore  be  written  in  the  form 


u 

V  " 


1 

4* 


if,  »te-) 


—  J  dxx  dyx 
c 


1 

4a 


z  dxx  dyx 


A2p  -  ;  === 

</[  (x— xi)2  —  (y-yi)2  -  z2]3 


(24) 


If  equation  (24)  is  transformed  to  the  original  apace  variables, 
the  relation  for  u  ia 
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Equation  (25)  is  valid  for  arbitrary  plan  forms  vith  known  load 
distributions.  Particular  examples  which  may  "be  worked  out  with 
relative  ease  are  the  lifting  surfaces  carrying  constant  load.  Once 
u  is  known  the  value  of  w  can  be  determined  fram#  the  integral 


bz  vi_oo 


and  from  w  the  ordinate  z  of  the  surface  as  a  function  of  x  and 
y  is  given  by 


where  Z.e.  denotes  the  leading  edge.  A  discussion  of  trapezoidal, 
rectangular,  and  triangular  plan  forms  with  constant  loading  is 
given  in  reference  2  although  the  method  of  derivation  is  different. 

Interest  in  constantly  loaded  wings  has  been  based  primarily 
on  the  fact  that  in  certain  cases  they  can  be  combined  to  produce 
surfaces  of  given  camber.  Thus,  a  superposition  of  trapezoidal  plan 
forms  of  variable  rake,  the  constant  loading  over  each  trapezoid 
being  a  function  of  its  rake  angle,  can  be  used  to  produce  a  flat 
plate  of  trapezoidal  or  rectangular  plan  form  at  an  arbitrary  angle 
of  attack.  In  this  case  the  loading  as  a  function  of  rake  angle  is 
determined  so  that  induced  vertical  velocity  is  kept  constant.  For 
problems  in  conical  flow  a  lifting  element  can  be  constructed  by 
subtracting  from  a  constantly  loaded  right  triangle  with  angle  of 
sweep  equal  to  6  the  constantly  loaded  right  triangle  with  sweep 
angle  equal  to  6  —  d5.  The  resultant  element  carries  a  constant 
load  and  has  a  sweep  angle  equal  to  6.  By  summing  these  elements 
it  is  possible  to  find  the  load  distribution  as  a  function  of  6 
such  that  certain  flat  lifting  surfaces  at  angles  of  attack  are 
formed.  In  reference  2,  triangular  wings  swept  back  of  the  Mach 
cone  were  studied  by  this  method  for  arbitrary  angles  of  yaw.  Brown 
(reference  7)  has  used  this  same  lifting  element  to  study  the  more 
restrictive  case  of  the  symmetrical  triangular  wing. 

A  brief  discussion  of  differences  existing  between  the  methods 
for  producing  the  swepi>-back  lifting  element  will  shed  some  light 
on  the  various  lines  of  attack.  The  approach  used  in  reference  2 
is  essentially  mathematical  in  that  a  particular  solution  of  the 
partial  differential  equation  is  used  in  conjunction  with  Green’s 
theorem  to  satisfy  the  boundary  conditions  of  the  problem.  The 
principal  criticism  of  such  a  method  is  that  the  physical 
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Interpretation  is  missing.  The  use  of  equation  (25) ,  however, 
removes  all  such  criticisms  for  precisely  as  in  the  case  of 
incompressible  flow  the  lifting  element  is  created  by  distributing 
doublets  over  the  wing.  In  Brown* s  solution  it  was  necessary  for 
him  to  deteimine  first  a  line  of  sources  by  means  of  an  integration 
along  the  line  and  then  to  form  the  doublet  line  by  differentiating 
along  the  normal  to  the  line.  The  order  of  differentiation  for 
incompressible  flow  is  immaterial,  since  the  limits  in  the  integral 
are  independent  of  the  position  of  the  point  PatX,  Y,  Z.  Super¬ 
sonic  flow  destroys  this  property  and  it  is  only  after  the  introduc¬ 
tion  of  the  concept  of  "finite  part"  that  the  derivative  of  an  integral 
may  be  written  as  the  integral  of  the  differential  coefficient  of  the 
Integrand*  Equation  (25)  thus  simplifies  the  analysis  and  at  the 
same  time  maintains  the  analogy  with  previous  work. 


Load  Distribution  for  Rolling  Wing 

The  usefulness  of  equation  (20)  is  not  at  all  restricted  to  a 
synthesis  of  previously  known  solutions.  As  an  example  of  its 
generality  consider  its  application  to  the  problem  of  the  rolling 
wing  with  leading  edge  swept  ahead  of  the  Mach  cone.  Figure  k 
shows  the  boundary  conditions  involved*  The  value  of  w  is 
specified  over  the  wing  and,  since  the  Mach  cone  is  behind  the 
leading  edge,  the  value  of  the  perturbation  velocities  u,  v,  and  w 
are  of  course  zero  ahead  of  the  leading  edge.  Assume  for  the 
moment  that  a  symmetrical  body  at  zero  angle  of  attack  is  considered. 
It  follows  that  if  =  vu  and  fl  j  =  Wj  then  equation  (20)  can 
be  written  in  the  form 


1 

w  =  — 
k 


(26) 


since,  for  reasons  of  symmetry,  the  normal  gradients  of  w  on  the 
two  surfaces  are  equal.  Using  now  the  fact  that  the  Mach  cone  is 
behind  the  leading  edge  then  the  pressure  over  the  upper  surface  is 
independent  of  the  shape  of  the  lower  surface  and  equation  (26) 
may  be  applied  directly  to  the  rolling  flat  plate  if  wu  is 
determined  from  the  given  induced  velocity  on  either  the  upper  or 
lower  surface.  This  method  of  approach,  of  course,  limits  the 
solution  to  cases  where  the  leading  edge  is  ahead  of  the  Mach  cone. 


If  the  rate 'of  roll  is  given  as  P  radians  per  second  then 
2wu  =  2PYi  and  equation  (26)  becomes 
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w 


w 


p 

X 


p*  YiZ  dYi  dXi 

t  (X— Xi) 2  -  p2  (Y— T1)2  -  P2  Z2  ]3/2 


(27) 


The  area  t  In  equation  (27)  is  that  contained  between  the  leading 
edge  and  the  trace  of  the  forecone  on  the  XY— plane.  Figure  5(a) 
shows  the  configuration  for  three  traces  corresponding  to  forecanes 
from  the  points  Pi,  P2,  and  P3  .  The  region  containing  the 
point  P2  da  distinguished  from  that  containing  Pi  and  P3  by 
the  fact  that  t  for  Pa  lies  ahead  of  the  Mach  cone  from  the  apex 
and,  furthermore ,  entirely  on  the  right  of  the  Xi  ails.  The 
regions  corresponding  to  Pi  and  P3  differ  in  the  fact  that  when 
integrating  from.  +00  to  z  to  find  u,  the  upper  limit  of  the 
integral  in  the  first  case  is  the  Mach  cone  X2  —  fi2Xs  —  p^2  =  0 
whereas  in  the  latter  case  the  upper  limit  is  the  leading-edge 

wedge  Z  =  -- — •  (fig.  5(t)). 
v  m2p2— 1 

The  solution  must  be  carried  out  separately  for  each  of  these 
regions  but  only  the  details  for  the  region  corresponding  to  Pi 
need  be  given  here  since  the  others  are  similar. 

It  follows  that  the  induced  velocity  u  at  the  point  P  is 
given  as  the  sum  of  the  triple  integrals. 


B  P  P 
u  =  ax  it 


/  dZQ  f  YidYi  P  y  - — — - -~T7Z 

'Wo>  J--£ 


p%dXi 


a_  p 
ax  * 


S 


A(-1,Z0) 


-p^/x^-p^ 


n*  P^odXi 

dZ0  /  YidYi  L  - - 


m 


=  y  a_p  rz 
I—*  ax  fl-4.  - - 

fc=— 1,1  P  *A2-p2Y2 


dZ0k 


[(X-Xi )  2-p2( Y-Yi )  ^-p2zo2  ]  3/2 


ZoYr  (X  +  k  dYi 


A(k,Zn)  t(Y-Yi)2+Z02]  J(T+%±) 2— p2(Y— Yi) 


(28) 
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where  A(k,  Z0)  is  the  value  of  Yi  determine d  hy  the  intersection 
of  the  forecone  with  the  leading  edge  on  the  right  and  left  sides 
with  k  equal  respectively  to  —1  and  +1.  See  fig.  5(a).)  Thus 


A(k,Z0) 


—  (k  jg  +  P3!)^  +  kfi  */(X  +  k  jj)2  +  Z02(s3  —  P2) 


After  reversing  the  order  of  integration  and  integrating  with 
respect  to  V  it  follows  that 


-  f 

bT  J. , 


A(k,Z) 


r(x+kj)  -  J ( X+k^) 2-p 2( Y-Yi )  2-p  *z  2 
YidYxZn  - . 1  j  .in  1 1 1  ■■r:,  u  a..a" - - 

L(X+k§0  +  J  (X+k|i)  2-*2( Y-Yx )  2-P^2 


Moreover,  since  the  integrand  is  zero  at  Yj.  =  A(k,Z)  the 
derivative  with  respect  to  X  can  be  taken  inside  the  integral  and 


I 


*i  dYi 


A(k,Z)  ,/(X  +  ~k)2  -  p2  (Y-Yi)2  -  p2  Z; 


Integrating  in  this  equation  and  combining  terms  it  follows  that 
induced  velocity  u  is  given  by  the  expression 


mp2!  +  X 
(m2P2— 1) 


P  2 
u  »  —  m 


arc  sin 


mp^  +  X 


P  2 


mB^Y  -  X 
(m2p2— l)3/2  L 


arc  sin 


J  (mX+Y)2  +  (1-rm2  P2)Z2 
mp2!  -  X 


( 


P^(mX—Y)a  +  (l-m2p2)Z2 


>f] 

y* . 


(29) 
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Setting  Z  =  0  in  equation  (29),  pressure  coefficient 
Ap  2u 

— -  * - is  given  by  the  expression 

q  v 


— 

q 


2P  f  mfi2Y  -  X  r  mp2!  -  X  1 

—  m2  i  - - —  arc  sin  . . — - 

i(7  L  (m2P2  -  1)  '  L  P(mX  -  T)  J 


mp2!  +  I 

(m2p2  -  1)  3/2 


euro 


r  mp2!  +  X  ‘ 

mp2! « 

L  p(mX  +  Y)  - 

+  (m2pa  -  1)  3/2J 

i}'(30) 


when  X>pY  , 

This  solution  holds  in  both  regions  containing  the  points  Pi 
and  Ps  .  However,  in  the  region  ahead  of  the  Mach  cone  hut  still 
on  the  wing  (region  corresponding  to  P2)  it  is  easy  to  show  that 


q 


2Pm2 

Y(m2p2  -  1)  a/2 


(mP2Y  -  X) 


(31) 


where  f-j-  <X<  PY.  Figure  6  shows  a  spanwise  plot  of  —  for 

m  q 

m  *  2  and  p  *  1. 


Equations  (30)  and  ( 31)  provide  sufficient  information  for  the 
calculation  of  the  stability  derivative  for  damping  in  roll,  Cj^. 

Integration  of  the  load  distribution  yields  the  result  that 


a(Pb/27) 


1 

3P 


Load  Distribution  for  Pitching  ¥ing 

Another  simple  application  of  equation  (20)  is  found  in  the 
solution  to  the  problem  of  the  pitching  wing.  Figure  7  shows  the 
boundary  condition  involved  which  is  that  the  vertical  induced 
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velocity  be  a  linear  function  of  Xx.  If  Q  is  the  rate  of  pitch  in 
radians  per  second*  then  w  =  on  the  wing.  Again  the  solution 

is  obtained  only  for  wings  which  have  leading  edges  swept  ahead  of 
the  Mach  cone.  (Although  solutions  can  he  obtained  for  leading  edges 
swept  behind  the  Mach  cone,  they  involve  integral  equations  and  do 
nothing  to  illustrate  further  the  direct  methods  of  this  report.) 

In  the  rolling  wing  case,  Q  was  set  equal  to  perturbation 
velocity  w  and  as  a  result  a  distribution  of  doublets  was  used 
in  equations  (27)  and  (28).  As  an  example  of  the  manner  in  which 
source— sink  distributions  may  be  used  for  the  same  type  of  problem, 
equation  (22)  will  be  applied  in  the  present  case.  Since  the  wing 
is  swept  ahead  of  the  foremost  Mach  cone,  induced  effects  on  the 
upper  and  lower  surface  are  independent  and 


Again  three  regions  containing  the  points  Pi,  P2*  and  P3  are 
distinguished  (fig.  5(a))  and  the  solutions  will  be  derived  only 
for  the  region  containing  Pi.  Integrating  first  with  respect  to 
Yi  and  then  differentiating  with  respect  to  X  yields 


B(k,Z) 


ks— 1,1 


Y-kmXj. 

arc  sin  - -  - —  ■  -  ■  -  dXj 

/(X-Xi)2  _ 

(33) 


where 


B(k,Z)  =  m 


X  -  p2!!™  -  pk  { [Xm-kY)2  +  Z2  (l-m2p2) 


1  -  m2p2 
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)-  (>  3  -1) 


q  '  2q 


Z^2r2)(m2p^.1)  _  ml(2-^2p2) 


mX(2-m2p2)-Y  /X-P^Ym' 

arc  sin 


VpY-£mX/ 


mX(2-tt2p2)+Y  /X+P2^ 

+ - arc  sin  - 


\  pY+BmX 


(34) 


Formula  (34)  Is  valid  for  the  regions  and  P3  of  figure  5(a). 

For  the  region  P2  the  solution  is: 


— )  ~r  (m2p2-l) 3/2  =  Y  -  2Xm  +  Xm3p2 


(35) 


Equations  (34)  and  (35)  provide  sufficient  information  for  the 
calculation  of  the  stability  derivative  for  the  damping  in  pitch  niq 

(£s&\  A/otA 

Values  of  — J  as  determined  from  these  equations  are 

plotted  in  figure  8  for  m  =  2  and  p  =  1. 
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APPENDIX  A 

Discussion  of  Conditions  on  Surface  X 

By  definition  X  is  the  surface  on  which  the  streamlines  of  the 
flow  first  experience  pressure  disturbances,  that  is,  the  surface 
along  which  the  stream  first  becomes  aware  of  the  existence  of  the 
wing.  Figures  2(a)  and  2(b)  were  introduced  to  show  the  nature  of 
the  configurations  involved  for  two  different  plan  forms.  It  is 
apparent  that  when  the  wing  is  swept  ahead  of  the  foremost  Mach  cone 
the  wedge-like  form  of  X  is  comparable  to  the  wedge  appearing  in 
purely  two-dimensional  problems  while  the  wing  swept  back  of  the  Mach 
cone  has  for  its  surface  a  conical  surface  and  thus  may  be  thought  of 
as  involving  a  purely  three-dimensional  problem. 
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In  order  to  determine  the  value  of  ®  on  X  it  is  sufficient 
to  impose  the  condition  that  the  tangential  component  of  velocity  is 
continuous  across  X  .  Such  a  condition  represents  no  essential 
restriction  since  it  is  an  immediate  consequence  of  continuity  of 
mass  flow  and  continuity  of  the  tangential  component  of  momentum 
across  the  surface.  As  a  result  of  this  condition,  however,  it 
follows 'that  the  tangential  component  of  the  perturbation  velocity 
is  zero  on  the  downstream  surface  of  X  since  it  is  obviously  zero 
on  the  upstream  surface.  Moreover,  velocity  being  equal  to  the 
gradient  of  the  velocity  potential  the  perturbation— velocity  potential 
must  be  equal  to  a  constant  on  X  .  But  an  arbitrary  constant  can  be 
added  or  subtracted  from  the  velocity  potential  so  that  with  no  loss 
of  generality  the  value  of  $  on  X  can  be  assumed  equal  to  zero 


and,  since  the  conormal  lies  on  the  surface  X 


^  is  also  zero. 


The  complete  analogue  to  equation  (10)  has  npw  been  developed  for 
&  =  $  so  that 


0(x,y,z) 


1_ 

2n 


rr\L.(^ 

J  Jj  Lrc  Mz 


dyi 


(Al) 

When  &  is  equal  to  one  of  the  perturbation— velocity 
components,  it  is  obvious  that  boundary  conditions  over  X  and  t 
cannot  be  considered  to  be  absolutely  arbitrary  since  it  Is 
necessary  to  include  the  added  restriction  that  the  resultant 
potential  ®  also  satisfies  the  equation 


sg$  a2$  a2$  _  Q 

bx2  by2  bz2 

Considering  first  the  case  where  the  wing  is  swept  behind  the  Mach 
cone,  it  follows  that 


ft  =  f  u(x1,y,z)dx1  =  P 

J , - 

vz2+y2 


y 

z2 


w(x,y,zi)dz! 
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» 


and,  after  evaluating  the  partial  derivatives  of  $  and  substituting 
in  the  given  differential  equation,  direct  calculation  leads  to  the 
conclusion  that  on  X  the  following  differential  equations  hold 


2j  +  2l  .  Ul 
07  oz 


=  0 


o  dvx  _  Stx 
2x  — =■  +  2z  — =■  +  Tx 

Bx  Bz 


=  0 


ai  2a  ♦  +  va  =  0 

Bx  By 


where  ux,  Vi,  and  Wx  are  the  valves  of  u,  v,  and  w  on  X. 

The  general  solutions  of  these  linear  partial  differential  equations 
can  he  written  as  follows 


ui 


7===  fi  (f ) » 

VyaTz2  ^z/ 


Wx 


It  has  been  stated,  however,  that  the  tangential  component  of  the 
total  perturbation—velocity  vector  vanishes  on  X  ,  or,  in 
analytical  terms 


Zux  +  mvi  +  nwx  -  0 

where  Z,  m,  n  are  direction  numbers  of  any  tangent  to  X  and 
therefore  satisfy  the  relation 

Zx  —  my  -  nz  =  0 

Substituting  the  known  expressions  for  ux*  Vi,  Wx,  it  follows  that 


or,  using  a  different  notation. 
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7 Fi  (f ) +  f  Fe  (S)  *  I  ^  <3 


Consider  now  the  special  case  when  Z  =  0  and  m 


■.  Under 


these  conditions 


bo  that 


(?)- 1  *4) 


Since  the  variables  ^  and  ~  are  separated,  the  solution  of  this 


equation  may  be  written 


?2  lx  )  ~  K  8114  F» 


where  K  is  a  constant.  Returning  now  to  the  case  where  l,  m,  and 
n  are  in  the  ratio  x:  y:  z,  direct  substitution  into  equation (A2) 
gives 


so  that 


5l!  I  +lI+ili-  0 

y  1  \zj  z  x  y 


71  *1)' 


_E3L(5.  +  JL 
1  Vy  z, 


and 
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This  equation  can,  however,  he  written  in  the  form 

£  +  Z  K  1 

y  z 


from  which  it  follows  that  K  =  0  and  Fi  =  F2  =  F3  =0.  All 
perturbation  velocity  components  are  seen  to  be  zero,  consequently 
equation  (Al)  is  valid  for  all  cases  in  which  the  wing  is  swept  back 
of  the  Mach  cone. 

A  discussion  of  conditions  on  the  surface  X  will  next  be  given 
for  the  case  where  the  leading  edge  of  the  wing  lies  along  the  y-exis 
(fig.  2(a))  and  where  Q  represents  u  or  w.  The  perturbation- 
velocity  potential  ®  may  be  given  by  the  relation 

u(xj.,y,z)  dxx 

J±z 

where  the  plus  and  minus  signs  in  the  limit  apply,  respectively,  to 
conditions  above  and  below  the  xy-plane.  Since  <b  must  satisfy 
the  basic  differential  equation,  an  added  restriction  is  Imposed  on 
u  and  as  a  result  of  this  condition  it  can  be  shown  that 


dz 


where  ui  is  the  value  of  u  on  either  the  lower  or  upper  surface 
of  the  wedge.  It  follows  that  the  values  of  u  on  the  two  surfaces 
are 


«i  *  fi(x,y)  and  ui  «  f2(x,y) 

and,  since  the  solution  is  independent  of  z,  and  x  is 
proportional  to  z  in  both  cases,  the  final  expressions  are 


ui  “  fi(y)  and  ui  -  f2(y) 


^  ®  (x>y»z)  is  defined  as  an  integral  involving  w,  the  same  type 
of  analysis  leads  to  the  conclusion  that  w  on  the  two  surfaces  can 
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also  be  expressed  as  functions  of  y  alone.  Perturbation  velooity 
v  will  not  be  considered  for  this  type  of  leading  edge  since  the 
inclusion  of  u  and  w  oovers  all  commonly  used  boundary  conditions. 

It  remains  to  substitute  the  results  Just  obtained  into  equa¬ 
tion  (19)  in  order  to  study  the  contribution  of  the  integrals  over 
X  .  Apparently  only  one  term  in  each  integrand  need  be  considered 
since  the  conormal  is  perpendicular  to  the  y— axis  and  the  gradient  of 
n  in  that  direction  vanishes.  As  a  preliminary  step  to  setting  up 
the  integrals  it  is  convenient  to  introduce  a  new  coordinate  system 
x" ,  y" j  z"  which  is  obtained  by  rotating  the  axial  system  about  the 
y-axis  so  that  the  x"—  and  zn-axes  lie  respectively  in  the  lower  and 
upper  wedge  planes  while  the  y"-axis  coincides  with  the  y-axis. 

The  transformation  of  variables  is 


x"  *  (Xi-Zi) 

•T? 


ZM  M  -  (Xi+Zi) 

S? 

When  ft  «  u,  the  last  two  integrals  in  equation  (19)  may  now  be 
written 


1 

S 


_i_ 

2* 


dx" 


Substituting  for  rc. 


this  expression  becomes 
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i_  r3*  filnUn 

ZtJj-J  ia-z2  V,z2-(y-y1)e-z2 


i_  p7+J  xs~zz  fg(yi)^yi 
2*  7yV  I^-Z2  <V/I^.(y_yi)2^2 


It  is  apparent  that  if  fi(y)  =  —  fa(y)  the  integrals  combine  to 
give  zero  so  that  equation  (Al)  may  again  be  used  in  all  calculations; 
moreover,  the  same  condition  applies  when  ft  =  v#  The  assumption 
that  fi(y)  -  —  f2(y)  is  equivalent  to  postulating  that  in  all  cases 
*i(y)  and  f2(y)  are  odd  functions  of  y.  In  application,  however, 
this  property  is  always  maintained. 

It  remains  finally  to  consider  the  case  when  the  leading  edge 
of  a  wing  is  swept  ahead  of  the  Mach  cone  and  when  ft  is  a 
perturbation-velocity  component.  As  a  means  of  avoiding  unnecessary 
complication  in  treating  the  problem  it  is  possible  to  substitute 
first  the  transformation  (reference  5) 


t  =  —  y  +  mx 
r\  *  —  x  +  my 

£  *  Vm^~l  z 

where  the  leading  edge  has  the  equation  y  =  mx,  z  =  0,  and 
m>l.  Direct  calculation  shows  that  the  basic  differential  equar- 
tion  and  the  Mach  cone  are  invariant  under  the  transformation  and 
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that  In  the  new  oblique  coordinate  system  the  leading  edge  lies 
along  the  n-axis  while  the  planes  of  the  wedge  became 

I  +  i  -  0 


Because  of  the  invariant  properties  of  the  transformation,  and  the 
fact  that  the  z  =  0  plane  is  fixed,  equation  (19)  is  applicable 
directly  to  the  boundary  value  problem  for  the  swept  wing  in  the 
new  coordinate  system.  The  treatment  of  the  integrals  over  X  can 
therefore  be  developed  algebraically  in  exactly  the  same  manner  that 
applied  to  the  previous  case,  hence  u  and  w  are  constant  along 
the  lines 


1  +  ^  =  0,  r\  =  conBt 

and,  again,  if  conditions  of  skew  symmetry  are  maintained  above  and 
below  the  z  =4=0  plane  the  integration  over  the  surfaces  X 
cancel.  Thus  equation  (20 )  is  seen  to  be  valid  for  0  equal  to 
perturbation-velocity  potential  or  perturbation  velocity  for  all 
types  of  striaght  leading-edge  configurations.  And  this  is  the 
complete  analogue  of  the  subsonic  theory. 
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n,  =  cos  a 
nz  -  cos  b 

n3  =  COS  C 


^  =  cos  (tt-gl) 
A  =  cos  A 
A  -  cos  c 


FIGURE  i-  GEOMETRIC  RELATIONS  BETWEEN  DIRECTION 

cosines  (n, ,n2,n^  of  normal  and  direction  cosines 
W>A\A)  OF  CONORMAL  TO  SURFACE  S. 
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fvgc/re  6  -  l/ar/of/on  of  pressure:  coef- 
f/c/en /  for  ro/Z/ng  Zr/ongu/ar  m’ng 
skvepZ  u/ieoc/  of  moc/i  cone . 


figure  <3  _  If? r/af /on  of  pressure  coef- 
f/c/enf  for  p/fc/i/ng  tr/ongu/or  W/hg 
swept  ufieuc/  of  muc/i  cone. 


